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Preface 


This  Final  Technical  Report  was  Prepared  by  the  Advanced  Engineering 
Technologies  Department,  Aircraft  Engine  Business  Group,  General  Electric 
Company,  Evendale,  Ohio  for  the  United  States  Air  Force  Systems  Command, 

Air  Force  Wright  Aeronautical  Laboratories  Wrlght-Patterson  Air  Force  Base, 
Ohio  under  Contract  F33615-81-C-2090.  The  work  was  performed  over  a  period 
of  two  and  one-half  years  starting  In  September  1961.  Luclen  L.  Debruge 
was  the  Air  Force  Project  Engineer  for  this  program. 

—  ^  The  objective  of  this  program  was  to  develop  and  codify  a  method  for 
predicting  the  spanwlse  redistribution  of  energy  and  loss  associated  with 
rotor/stator  wakes  In  an  axial-flow  compressor.  The  mechanisms  considered 
were? (1)  spanwlse  transport  of  wake  fluid  due  to  Imbalance  of  the  radial 
pressure  gradient^ and  (2)  accumulatlon/dllutlon  of  wake  fluid  at  the  Inner 
and  outer  casing.  The  work  consisted  of  /developing  computer  modules  that 
perform  the  above  defined  objective. 

For  the  General  Electric  Company  Mr.  J.S.  Keith  was  the  Technical 
Program  Manager  for  this  program.  Dr.  C.E.  Whitfield  was  the  principal 
Investigator.  Mr.  A.J.  Bllhardt  was  the  overall  Program  Manager. 
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Equation  (18)  can  be  simplified  further.  Writing 


dU].  -►  Duj.  Duj-  ^  . 

— -  +  w.Vu  *  —  =  W  ■=: —  (x  IS  streamwise) 

3t  r  Dt  Dx 


vq  *  Vq  *  2Qr  ■  Cg  +  cg 

where  Cg  is  the  absolute  tangential  velocity  in  the  freestream,  cg  is  the 
absolute  tangential  velocity  in  the  wake  and 

ue  “  *'0  ■ 

■  (wx  -  W)  sin  B 

■  D  W  F  sin  B  (from  the  previous  section). 

Neglecting  the  fourth  and  last  terms  in  the  left-hand  side  of  Equation  (18), 
in  accordance  with  Equation  (1)  of  Section  3.1.1,  and  introducing  the  concept 
of  eddy  viscosity,  as  was  done  in  the  previous  section,  gives  the  following 
"simplified  radial  momentum  equation"  for  the  wake  fluid: 


Ce  + 

r 


D  W  F  sin  6  +  e 


(19) 


The  value  of  cg  in  Equation  (19)  can  be  calculated  from  the  expression 


c„  =  W(1  -  DF)  sin  6  +  (^r 

y 


The  results  given  in  Reference  12  were  obtained  using  the  approximation 

Cg  =  C  g 


which  is  felt  to  be  reasonable  in  the  far  wake  except  when  the  absolute 
velocity  is  near  axial. 

Assume  that  the  radial  velocity  profile  of  the  wake  can  be  represented 
by  the  same  Gaussian  function  as  the  streamwise  profile 


W 


"r  -  ^r 
W 


G  F  -  G  e 


-wn^ 
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Here  w  is  the  velocity  of  any  wake  particle  and  Fy  represents  the  viscous 
forces  in  the  wake.  The  wake  incremental  velocity,  u,  is  defined  by: 


w 


u  +  W 


%»here  W  is  the  freestream  velocity,  and  will  denote  the  freestream 
fluid  density. 

Substituting  for  ti  in  Equation  (16)  gives 


+  w'Vu  +  u*7W  +  2n  X 


“  ^  (pt '  t)''"  - 


-{f- 


•VW  +  2^  X  W  -  + 


■  ° 


(17) 


The  right-hand  side  of  the  above  equation  is  the  equation  of  motion  for  the 
freestream  fluid  which  is  assumed  to  be  satisfied. 

Using  cylindrical  polar  coordinates  and  making  use  of  the  vector  rela¬ 
tionships 


3(i9)  .  .  1 

36  r 


and 


3(ir) 

36 


gives,  for  the  radial  component  of  Equation  (17) 


3u 

“3t 


_  +  .  __  +  -  —  -  2nu^  +  j 


(18) 
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Substituting  this  into  Equation  (12)  gives  the  results 


Therefore,  if  the  eddy  viscosity,  e,  and  the  freestream  velocity,  W,  are  con¬ 
stant  with  distance  dotmstream,  the  wake  thickness  and  velocity  deficit  will 
vary  as 


6  *  (x  -  Xq)' 


(15a) 


D  «  (x  - 


(15b) 


The  choice  of  Xq  in  the  results  presented  in  References  11  and  12  was  deter¬ 
mined  empirically  from  experimental  data  as  a  point  12Z  of  chord  upstream  of 
the  trailing  edge. 


3.2.2  Radial  Wake  Decay 

A  similar  approach  is  used  to  calculate  the  decay  of  the  peak  radial 
velocity  in  the  wake  with  distance  downstream.  The  equation  of  motion  for 
flow  in  rotating  coordinates  is 


+  (w*V)w  +  2Qxw  -  r  +  “  Fy 


Knowing  the  form  factor  at  the  trailing  edge,  Hte>  using  the  above,  the 
maximum  streamwise  deficit  at  the  trailing  edge,  D^g,  is  given  by 


TE 


(10) 


3.2  WAKE  DEVELOPMENT 

3.2.1  Streamwise  Wake  Decay 

In  order  to  follow  the  path  of  the  wake  fluid,  it  is  necessary  to  predict 
the  decay  of  the  wake  velocity  defect.  This  is  accomplished  by  considering 
the  streamwise  momentum  equation  for  the  wake  which  can  be  written 


w 


X 


+  w 


3w 
_ X 

y  3y 


+  i  ^ 

p  dx 


t 


3y2 


(11) 


In  the  "far  wake"  (greater  than  some  fractional  blade  chord  distance  down¬ 
stream),  Wjj  s  Also,  from  the  freestream  momentum  equation 


i  ^ 

p  dx 


Neglecting  the  second  term  in  Equation  (11)  and  assuming  that  the  pressure 
and  density  in  the  wake  are  equal  to  the  freestream  values.  Equation  (11) 
becomes 


W 


3w 


_  w  iw 

3x  3x 


(12) 


where  e  is  an  eddy  viscosity. 

In  Section  3.1.2  the  wake  profile  was  described  by  a  universal  function: 


D  F(n) 


10 


CD  1*0 


and  its  momentum  thickness,  likewise 


00 

-/  ^  f  ■  oy-An 

—00  _  J 


Integrals  of  the  Gaussian  profile  function  are  found  to  be: 


Gii  -  If  dn  -  1 


Gi2  -y  F^on  -  J- 
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■/ 


F-^dn 


Thus,  the  wake  form  factor  is 
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1  -  D  ^12 

G,, 


1  -  -  D 

/T 
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Formulas  for  calculating  these  velocity  ratios  may  be  found  in  Appendix  1  of 
Reference  7. 

Given  the  value  of  Hxg,  Equations  (4)  and  (3)  are  substituted  into  Equa¬ 
tion  (2)  which  is  integrated  from  midchord  to  the  trailing  edge  of  the  blade. 
It  is  assumed  that  at  midchord  the  blade  boundary  layer  is  collateral.  It  is 
further  assumed  that  the  result  of  the  integration  of  Equation  (2)  is  a  radial 
velocity  increment  at  the  trailing  edge,  which  is  to  be  added  to  the  free- 
stream  radial  velocity.  Note  that  the  value  0.05  was  selected  so  Equation  (2) 
would  yield  the  peak  wake  radial  velocity  observed  experimentally. 


3.1.2  Streamwise  Wake  Velocity  Defect 

The  streamwise  wake  velocity  defect  is  based  on  the  assumption  that  all 
wake  profiles  can  be  represented  by  a  single  function 


W-w 

*  -  D  F(n) 


(6) 


where 


W  *  freestream  relative  velocity 
w^  *  streamwise  velocity  in  the  wake 

W-Wjj  Up 

D  ■  maximum  velocity  defect,  — 


F(ri)  =  Gaussian  profile  function  = 

"  ■  f 

5  ■  wake  thickness  parameter,  approximately  equal  to  the  wake  half  width 

x,  y  *  streamwise  and  cross-stream  coordinates. 

The  displacement  thickness  of  the  wake  is 


F  (n)  dn 


(7) 


8 


1 

k  -  (0.05)“  - 
w 


H-1 

(0.05)  ^ 


(3) 


In  order  to  integrate  Equation  (2)  along  the  blade  chord  we  must  find  an 
expression  for  H,  the  form  factor  of  the  streamvise  boundary  layer.  Assume 
that  boundary  layer  growth  is  most  pronounced  over  the  trailing  half  of  the 
blade.  Let  the  streamwise  boundary  layer  be  represented  by 


1 


at  the  midchord  of  the  blade.  Then  H  ~  >  1.25  at  this  point  and,  if  the 

chordwise  distribution  of  H  is  given  by  a  parabola,  we  have 


H(z) 


H  +  4 


«o> 


(z  - 


^o> 


/(zte  ' 


(4) 


where  Zq  is  at  midchord  and  Hxe  has  yet  to  be  determined. 

The  loss  correlations  of  Koch  and  Smith  (Reference  7)  show  the  variation 
of  trailing  edge  form  factor  and  momentum  thickness  with  equivalent  diffusion 
ratio.  In  this  work  we  are  interested  in  the  value  of  trailing  edge  form 
factor  resulting  from  the  suction  surface  only;  the  Koch  and  Smith  curves  have 
been  modified  to  reflect  this.  Since  only  leading  and  trailing  edge  informa¬ 
tion  is  assumed  known,  the  calculation  of  equivalent  diffusion  ratio  is 
related  to  the  following  three  velocity  ratios: 


V 

max 


'^oTE 


(5) 


where 

Vp  *■  mean  passage  velocity  in  the  throat  region, 
^max  ■  maximum  velocity  on  the  suction  surface, 

Vq  «  free  stream  velocity  before  wake  mixing 
Vi  *  inlet  velocity. 
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3.0  ANALYSIS  AND  COMPUTATIONAL  TECHNIQUES 


3.1  TRAILING  EDGE  HAKE  DEFINITION 


3.1.1  Wake  Peak  Radial  Velocity 


The  wake  radial  velocity  calculation  is  based  on  Equation  (32)  of  Adkins 
and  Smith  (Reference  1). 


where 


l^tan^  $2 

+  2 


tan  02 


fir 

Wz 


(1 


(1) 


Wz 


If  it  is  assumed  that  DW^.  ■  0  (no  free-stream  radial  acceleration),  Equation 
(1)  becomes 


-Dz  tan  $2 
r 


[tan  02 


(2) 


The  boundary  layer  streamwise  velocity  profile  can  be  approximated  by  a  power 
law  distribution; 


1 

n 


which  yields 

„  5  *  n+2 

”  “  T  “  “ 


It  is  reasonable  to  assume  that  conditions  in  the  boundary  layer  are 
typified  by  those  occurring  at  some  point  relatively  deep  in  it,  say  at 
y/6  *  0.05.  It  can  be  seen  (for  example,  Johnston  in  Reference  14)  that  this 
is  sufficiently  far  away  from  the  wall  for  inviscid  theory  to  be  used  in  esti¬ 
mating  the  radial  crossflow.  Thus, 


2.3  CONTENTS 


This  report  describes  the  techniques  used  to  predict  the  wake  peak  radial 
velocity,  the  peak  streamwise  velocity  decrement,  and  the  momentum  thickness 
and  displacement  thickness  at  the  trailing  edge  of  a  blade.  Also  defined  are 
the  "mass-energy-average"  velocity  decrements  of  the  wake  and  the  decay  of 
these  various  velocity  decrements  with  distance  downstream.  The  radial  veloc¬ 
ity  of  the  wake  fluid  will  also  be  affected  by  its  proximity  to  the  annulus 
walls.  Taking  these  factors  Into  account,  the  tracking  of  wakes  Is  described. 
Finally,  the  insertion  of  the  wake  source  terms  into  the  mixing  equation  and 
the  subsequent  effects  on  the  circumferential  average  flow  solution  are  shown. 

The  results  obtained  during  the  course  of  this  investigation  show  the 
agreement  between  theory  and  experimental  data  obtained  for  peak  radial  and 
streamwise  velocity  increments/decrements  in  the  wakes  of  four  machines, 
together  with  differences  in  peak  radial  velocities  calculated  by  the  Adkins 
and  Smith  method  and  by  the  method  described  here.  Comparisons  of  circumfer¬ 
ential  average  flow  solutions  obtained  with  the  Adkins  and  Smith  scheme,  and 
with  the  new  scheme  showing  both  the  effect  of  the  new  wake  definition  on  the 
mixing  results  and  the  effect  of  the  inclusion  of  wake  fluid  migration  in 
addition  to  this,  have  been  made  for  the  Air  Force  1500  ft/s.  Transonic,  High- 
Through-Flow,  Single-Stage  Axial-Flow  Compressor. 

Not  presented  in  this  final  report  is  the  effect  of  the  wake  impinging  on 
the  pressure  side  of  downstream  blade  rows  in  altering  the  radial  velocity  of 
the  wake.  This  effect  was  studied  and  then  omitted  because  it  was  found 
inconsequential . 


that  at  the  midchord  of  the  blade  the  boundary  layer  is  collateral  and 
the  form  factor  grows  parabolical ly  from  midchord  to  the  value  predicted  by 
the  modified  Koch  and  Smith  curves  at  the  trailing  edge  of  the  blade. 


These  assumptions  enable  the  Adkins  and  Smith  equation  for  the  peak 
radial  velocity  of  a  representative  particle  in  the  blade  suction  surface 
boundary  layer  to  be  integrated  from  the  midchord  to  the  trailing  edge  of  a 
blade.  When  combined  with  the  wake  radial  velocity  increment  decay  laws  form¬ 
ulated  in  the  course  of  this  work  and  described  here,  encouraging  agreement  is 
found  between  theory  and  experimental  data  for  most  test  cases  studied. 


2.2  MIXING  CALCULATION 


The  objective  behind  a  better  calculation  of  the  peak  radial  velocity  in 
the  wake  is  to  include  the  effects  of  wake  migration  in  the  movement  of  fluid 
properties  across  streamlines  -  the  mixing  process.  In  the  original  work  of 
Adkins  and  Smith,  wake  radial  velocity  and  thickness  were  used  to  modify  the 
spanwise  secondary  flow  profile  resulting  from  the  various  Tref ftz-plane  solu¬ 
tions  of  Poisson's  equation  for  the  secondary  flow  stream  function. 

In  the  work  described  here,  the  definition  of  the  wake  radial  velocity 
and  width  have  been  improved.  In  addition,  the  effect  of  wake  fluid  property 
migration,  as  opposed  to  mixing  of  circumferential  average  properties,  has  been 
modeled . 

Consider  a  wake  shed  from  a  rotor  blade.  The  fluid  in  the  wake  will  have 
a  lower  relative  streamwise  velocity  than  the  surrounding  fluid.  Because  of 
this,  the  wake  fluid  will  have  a  higher  absolute  swirl  velocity  and  hence, 
assuming  the  total  rothalpy  in  the  wake  is  equal  to  that  in  the  free  stream, 
a  higher  total  enthalpy  and  angular  momentum.  The  entropy  of  the  wake  fluid 
will  also  be  higher  than  that  of  the  surrounding  fluid.  As  the  wake  fluid 
progresses  downstream  it  will  carry  these  property  increments  with  it;  thus, 
as  the  fluid  migrates  radially,  there  will  be  a  redistribution  of  properties. 

In  order  to  calculate  the  effect  of  wake  fluid  migration  on  the  radial 
distribution  of  circumferential  average  properties,  it  is  necessary  to  define 
the  "property  carrying  capacity"  of  the  wake,  the  streamwise  and  radial  veloc¬ 
ities  at  which  the  properties  will  be  carried,  and  the  radial  location  of  the 
wake  fluid  at  the  various  calculation  stations  downstream  of  the  shedding  row. 
To  include  this  property  migration  in  the  circumferential  average  flow  solu¬ 
tion,  it  is  necessary  to  modify  the  mixing  equation  of  Adkins  and  Smith  to 
include  these  wake  "source"  terms  and  use  this  new  form  of  the  equation  in  the 
through-flow  computer  analysis. 
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The  calculations  of  Adkins  and  Smith  demonstrate  quite  dramatically  the 
effects  of  enabling  and  disabling  the  mixing  process.  Agreement  obtained 
between  calculations  and  test  data  with  mixing  enabled  showed  clearly  the  need 
for  taking  this  process  into  account  when  calculating  the  circumferential 
average  flow  solution.  However,  at  the  outset  of  the  work  described  here  it 
was  felt  that  the  blade  boundary  layer/wake  centrifugation  model  used  by 
Adkins  and  Smith  was  considerably  less  sophisticated  than  their  other  work. 
First,  it  was  felt  that  the  trailing  edge  estimates  of  wake  width  and  the 
spanwise  component  of  velocity  were  unnecessarily  crude.  Second,  their  model¬ 
ing  of  the  "mixing"  process  through  the  use  of  the  homogeneous  diffusion  equa¬ 
tion  precluded  the  possibility  of  predicting  the  tendency  for  the  high  tem¬ 
perature  fluid  in  the  rotor  wake  to  migrate  outward,  accumulate  near  the  cas¬ 
ing,  and  increase  the  circumferential  average  temperature  there. 

Thus,  the  work  described  here  and  in  References  8  through  13  falls  into 
two  categories:  improvement  of  the  blade  boundary  layer/wake  centrifugation 
model  of  Adkins  and  Smith  to  bring  it  to  the  same  level  of  sophistication  as 
the  other  four  secondary  flow  models,  and  reformulation  of  their  mixing  model 
in  order  to  present  a  more  realistic  picture  of  the  flow  properties  as  they 
are  convected  downstream. 


2.1  WAKE  RADIAL  VELOCITY  DEFINITION 

The  method  used  by  Adkins  and  Smith  to  determine  the  maximum  radial 
velocity  in  the  wake  shed  from  a  rotating  or  stationary  blade  in  an  axial-flow 
turbomachine  is  believed  to  be  deficient  in  that  no  account  was  taken  of  the 
effects  of  blade  loading.  The  approach  taken  here  is  Intended  to  rectify  this 
situation. 

In  order  to  derive  an  expression  for  spanwise  velocity  in  the  boundary 
layer  on  an  axial-flow  turbomachine  blade,  Adkins  and  Smith  considered  the 
spanwise  acceleration  of  a  representative  small  mass  of  fluid  in  the  blade 
boundary  layer  at  a  representative  point  along  the  chord.  In  the  absence  of 
viscous  stresses,  this  acceleration  was  assumed  to  act  over  the  time  taken  by 
the  particle  to  travel  a  representative  distance,  resulting  in  a  spanwise 
velocity.  Comparison  with  available  test  data  led  to  the  selection  of  the 
value  of  a  constant  in  the  model. 

For  the  present  work,  the  equation  derived  by  Adkins  and  Smith  has  been 
retained.  Application  of  this  equation  requires  an  integration  along  the  chord 
of  the  blade  and  it  is  the  underlying  assumptions  behind  this  integration  that 
have  been  changed.  In  particular,  blade  loading  effects  have  been  included  by 
means  of  the  loss  correlations  of  Koch  and  Smith  (Reference  7)  which  have 
been  modified  for  this  purpose  to  reflect  boundary  layer  growth  on  the  blade 
suction  surface  only.  It  has  been  assumed  that  the  streamwlse  velocity  pro¬ 
file  of  the  boundary  layer  on  the  blade  suction  surface  can  be  represented 
by  a  power  law  distribution,  and  also  that  conditions  in  the  boundary  layer 
are  typified  by  those  occurring  at  a  point  relatively  deep  in  it  whose  location 
can  be  chosen  to  approximate  the  peak  radial  velocity.  It  is  assumed  further 


2.0  INTRODUCTION 


In  the  study  of  the  complex  flow  that  occurs  within  the  compressor  and 
turbine  of  today's  axial  flow  turbomachines  it  is  usual  for  simplifying 
assumptions  to  be  made.  For  instance,  virtually  all  analytical  methods  used 
at  present  in  the  design  of  compressor  and  turbine  blade  rows  contain  the 
assumption  that  the  flow  remains  on  axisymmetric  stream  surfaces  as  it  passes 
through  the  blade  row;  an  assumption  that,  as  designers  are  aware,  is  not 
really  valid. 

The  causes  of  blade-to-blade  stream  surface  distortion  are  many  and  vari¬ 
ous.  For  instance,  if  the  incoming  flow  to  a  blade  row  contains  vorticity  or 
if  the  blade  circulation  varies  along  the  span,  secondary  flows  are  generated 
that  cause  blade-to-blade  (or  SI)  stream  surface  twists.  Likewise,  if  a  sweep 
condition  is  present,  the  tendency  of  the  flow  to  maintain  its  spanwise  veloc¬ 
ity  component  while  other  velocity  components  are  changed  by  the  turning 
action  of  the  blades  will  lead  to  SI  surface  distortions.  Another  physical 
reason  why  some  portion  of  the  fluid  passing  through  a  blade  row  may  depart 
substantially  from  otherwise  axisymmetric  or  slightly  twisted  SI  surfaces 
involves  three  dimensional  flow  in  the  blade  boundary  layers  and  wakes.  It  is 
this  last  phenomenon  with  which  the  present  work  is  concerned. 

In  a  recent  paper,  Adkins  and  Smith  (Reference  1)  describe  a  method 
whereby  the  secondary  flows  generated  within  the  blade  passages  of  an  axial- 
flow  turbomachine  can  be  quantified  and  used  in  an  analysis  of  the  meridional 
flow  field.  In  the  course  of  this  work,  they  considered  that  the  secondary 
flow  field  was  the  resultant  of  five  effects:  (1)  mainstream  nonfree-vortex 
flow,  (2)  end  wall  boundary  layers,  (3)  blade  end  clearances,  (4)  blade  end 
shrouding,  and  (5)  blade  boundary  layer  and  wake  centrifugation.  The  tech¬ 
niques  used  were  based  on  earlier  work  of  Smith  (References  2  through  5)  and 
Lakshimlnarayana  and  Horlock  (Reference  6).  In  particular,  models  were  devel¬ 
oped  to  predict : 

a.  The  overturning/underturning  swirl  angle  deviation  that  must  be 
added  to  the  primary  (two-dimensional  blade-to-blade)  flow  to  arrive 
at  the  full  three-dimensional  flow  swirl  angles 

b.  Spanwise  redistribution  of  (circumferential  average)  total  tempera¬ 
ture  and  total  pressure 

c.  End  wall  losses  that  may  be  applied  on  a  blade-row  by  blade-row 
basis  and  used  with  the  profile  loss  model  of  Koch  and  Smith 
(Reference  7) . 

The  present  effort  is  directed  toward  improvement  of  Part  (b)  of  the  Adkins 
and  Smith  (Reference  1)  work.  Parts  (a)  and  (c)  are  l_ft  essentially  intact. 


1.0  SUMMARY 


This  report  describes  work  accomplished  under  Air  Force  Contract  FSSblS- 
81-C-2090  between  September  1981  and  March  1984.  During  this  time,  formulae 
have  been  derived  and  encoded  for  definition  of  the  peak  radial  and  stream- 
wise  velocity  decrements /increments  in  a  wake  at  a  blade  trailing  edge  and 
their  decay  with  downstream  distance.  Comparisons  with  data  show  good  agree¬ 
ment  in  most  cases.  In  addition,  the  thickness,  displacement  thickness,  and 
momentum  thickness  of  the  wake  have  been  calculated  and  used  in  the  prediction 
of  changes  in  the  circumferential  average  flow  solution  caused  by  the  migra¬ 
tion  of  wake  fluid  across  streamlines.  Wake  fluid  has  been  tracked  from  its 
shedding  blade  row  through  the  downstream  row  with  the  effects  of  fluid  buildup 
on  hub  and  casing  walls  taken  into  account.  Final  results  compare  circumfer¬ 
ential  average  flow  results  obtained  using  the  new  wake  definition  siodel  with 
earlier  analysis. 


.  SYMBOLS  (Concluded) 


Definition 

Free  stream 
Wake  aggregate 
Meridional 

Peak  or  centerline  value  within  the  wake 

Radial 

Spanwlse 

Trailing  edge 

Tangential 

Wake 

Streamwise 

Axial 

Origin 

Inlet  to  cascade 


Outlet  from  cascade 


where  G  is  the  maximum  value  of 


Substitution  in  Equation  (19)  leads  to  the  result 


dG  ^  2D  C9  p.  G 

dx  rW  2(x  -  Xq) 


if  W  is  constant  with  distance  downstream  and  Ce  ■  ce. 
equation  is 


In  integral  form  this 


G2  *  G 


^  A  L 


2Cq  sin  6 
"rW 


D  + 


2(x  -  Xq) 


dx 


(20) 


where  1  and  2  are  adjacent  points  in  the  streamwise  direction.  Equation  (20) 
is  integrated  numerically,  no  analytic  result  having  been  found. 


3.3  WAKE  PARTICLE  TRACKING 


3.3.1  Calculation  of  Wake  Fluid  Trajectories 

Consider  a  blade  row  in  an  axial-flow  turbomachine.  There  is  flow  rela¬ 
tive  to  this  row;  thus,  a  wake  is  shed  from  the  trailing  edge  of  each  blade  in 
the  row.  At  each  radial  streamline  position,  the  initial  deficits  in  the  wake 
in  both  streamwise  and  radial  directions  can  be  calculated  by  the  formulae  of 
Section  3.1.  As  the  wake  passes  downstream  we  can  calculate  the  decay  of  these 
deficits  by  the  application  of  formulae  from  Section  3.2.  Primarily  we  are 
concerned  with  the  radial  movement  of  the  wake  fluid  relative  to  the  circum¬ 
ferential  average  flow  streamlines.  This  is  accomplished  by  integrating 


•'TE 


2iTr  X  pC„  u 

m  rg 


cos  Dt, 


where  il^g  is  the  stream  function  value  of  the  wake  aggregate,  Ufg  is  the  radial 
velocity  increment  of  the  wake  fluid  (uj-g  *  w^g  -  W^),  Wtf  is  tne  total  mass 
flow  rate  and  X,  P,  Cg,,  and  9  are  the  circumferential  average  through-flow 
values  of  blockage,  density,  meridional  velocity,  and  meridional  angle. 


^  ^  pi  IP  1^1 


The  radial  velocity  increment,  ure,  may  take  into  account  not  only  the 
radial  velocity  in  the  wake  but  any  effects  of  secondary  flow  in  downstream 
blade  rows  and  any  modification  to  the  wake  radial  velocity  caused  by  contact 
with  downstream  blade  surfaces.  Such  effects  were  considered  in  some  detail 
in  the  early  part  of  this  study,  and  the  results  were  reported  in  References 
9,  10,  and  11.  However,  because  these  effects  are  small  (in  terms  of  the 
overall  mixing  process),  it  was  felt  that  the  additional  complexity  needed  to 
include  the  impingement  and  secondary  flow  modification  to  radial  wake  veloc¬ 
ities  in  the  product  code  was  not  justified  and,  hence,  has  been  omitted. 


3.4  TRANSPORT  OF  WAKE  PROPERTIES 
3.4.1  General 

Let  P  be  any  property.  Further,  let  P  be  the  circumferential  average  of 
P  and  AP  the  (mass-averaged)  increment  of  P  associated  with  a  spanwise  migrat¬ 
ing  "hot"  wake.  Specifically, 


AP 


dy 


Consider  the  spanwise  transport  of  P  when: 

1.  There  exists  a  blade-to-blade  variation  of  the  spanwise  velocity 
and  a  radial  variation  of  P,  but  AP  ■  0  everywhere,  or 

2.  There  is  a  radial  migration  of  the  "hot"  wake,  for  which  AP  is 
nonzero,  but  P(r)  is  radially  constant. 

These  two  cases  are  schematically  depicted  In  Figure  1. 

Conceptually,  both  of  these  cases  can  be  treated  by  the  inviscid  convec¬ 
tion  equations.  First,  identify  a  large  number  of  fluid  elements  at,  say,  the 
trailing  edge  plane,  and  then  use  a  (presumed)  knowledge  of  the  three- 
dimensional  flow  field  velocities  (especially  the  spanwise  velocity  components) 
to  compute  the  relative  position  of  these  fluid  elements  at  a  downstream  plane. 
From  the  conserved  convected  properties,  stagnation  enthalpy  (rothalpy)  or 
entropy,  of  each  fluid  element  a  determination  of  the  downstream  circumferen¬ 
tial  average  property  is  possible. 
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In  practice,  Adkins  and  Smith,  in  Section  4.1  of  Reference  I,  found  that 
the  spanwise  transport  of  P  (Case  1)  is  modeled  simply  and  accurately  by  the 
diffusion  equation*: 


where 


=  £sl  r 
aCm  J  s 

/s 


Unfortunately,  this  equation  does  not  model  a  convection  process  in  which  there 
exists  a  correlated  cross-passage  variation  of  property  P  and  radial  velocity 
w^;  that  is,  when  the  cross-passage  P  and  w^  profiles  have  roughly  the  same 
shape  (Case  2). 

A  method  for  treating  this  second  type  of  spanwise  transport  is  the  sub¬ 
ject  of  the  present  section.  The  method  for  combining  the  two  cases  is  subse¬ 
quently  described  in  Section  3.5. 

Before  proceeding,  a  point  of  clarification:  As  illustrated  in  Figure  2, 
one  can  define  two  values  of  AP 


AP  = 
w 


NgJ  max  (0,  P-P)  pw^  dy 
2iTrA  p  Cm 


Ng  I  min  (0,  P-P)  pw^^  dy 


2TrrX  p  C 


^Equation  (21)  was  derived  by  considering  only  spanwise  secondary  flow.  The 
equation  is  not  strictly  valid  in  a  region  with  a  large  cross-passage  compo¬ 
nent  of  secondary  velocity,  as  occurs  near  the  end  wall. 
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where  the  integrals  are  carried  out  over  one  blade  spacing  and 


aax  (A,  B)  s  larger  of  A  and  B 
min  (A,  B)  =  smaller  of  A  and  B 
From  the  above  it  follows  that: 

^w  ♦  -  0 

So,  APfg,  the  property  (or  energy)  increment  associated  with  the  freestream 
fluid,  is  the  same  magnitude  as  that  of  the  wake  fluid  ( AP^) .  However,  the 
radial  velocity  of  the  freestream  fluid,  opposite  in  sign  from  that  of  the  wake 
fluid  to  maintain  continuity,  is  much  smaller  in  magnitude.  Hence,  the  spanwlse 
displacement  of  energy  APfg  is  smaller  than  that  of  AP^  and  is  neglected. 

Furthermore,  it  is  convenient  to  use  the  freestream,  rather  than  the 
average  value,  as  the  reference  level.  Thus,  in  the  following  sections,  the 
freestream  value,  P^ <  will  replace  P  in  the  equations  that  define  AP. 


3.4.2  Passage-Mass-Weighted  Wake  Property  Increments 

As  described  in  Section  3.3,  particles  of  wake  fluid  are  generated  at  the 
trailing  edge  of  a  blade  row  and,  as  described  in  the  same  section,  are  tracked 
through  the  next  downstream  row.  These  particles  of  fluid  are  carrying  energy 
in  the  form  of  total  enthalpy  or  total  rothalpy,  and  they  also  have  entropy 
associated  with  them.  Therefore,  they  do  not  move  downstream  at  the  peak 
streamwise  velocity  decrement  calculated  earlier,  but  rather  at  a  mass-energy- 
decrement-weighted  average  velocity  obtained  as  follows.  Consider  an  aggre¬ 
gate  of  wake  fluid:  this  aggregate  has  mass 

♦  6 

Mg  -  At  ha  j  pwx  dy  (22) 

-A 

where 

ha  ~  lamina  thickness 

At  >  arbitrary  time  interval 


Assume  the  aggregate  is  carrying  energy  in  the  form  of  total  enthalpy  and  is 
part  of  a  wake  shed  from  a  rotor  blade.  Assume  also  that  the  total  rothalpy 
across  the  wake  is  constant  and  equal  to  the  freestream  value. 


Thus  in  the  wake 


H  -  I  +  U  ce^  -  I  +  U(U  +  we) 

H  -  I  +  U2  +  sin  B 

and  for  the  freestream 

Hg  «  I  +  U2  +  UW  sin  B 

where  B  is  generally  negative  for  a  rotor.  The  local  energy  increment  of  the 
wake  fluid  is,  thus, 

H  -  Hg  -  U(wx  “  W)  sin  B  (23) 

and  the  integrated  mass-energy-decrement  is 

(5 

Mg  6Hg  -  At  h^  J  pw^  dy 

6 

-  -At  h  p  UW^  sin  6  /  A  \\  ^ 

-6  e 

M  6H  -  -At  h  p  UW^  sin  3  0  (24) 

g  g  m  ^e 


where  6  is  the  momentum  thickness  in  the  wake. 

The  total  mass  of  fluid  for  a  given  lamina  thickness,  hQ,  is  given  as 


At  2TrrA 


N 


B 


m  ^  m 


(25) 


where  P  and  Cg,  are  circumferentially  mass-averaged  values  used  in  the  through- 
flow  analysis  program  and  Ng  is  the  number  of  blades.  From  this  it  follows 
that  the  passage-mass-weighted-energy  increment,  for  any  spanwise  position,  is: 


AH 


M  6H 


8  g 


UW^  sin  B 

X  C  2Trr 

m 


(26) 
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The  small  difference  between  the  freestream  density,  Pe,  and  the  average  dens¬ 
ity,  p,  has  been  neglected. 

In  addition  to  the  increment  of  energy,  the  wake  fluid  also  convects  an 
increment  of  entropy  (loss).  By  definition, 

S  -  -  C  In  (|  )  -R  In  (P  ) 

e  P  Te  p 

e 


where  S  is  the  entropy  at  any  point  in  the  wake  and  Sg  is  the  entropy  in  the 
freestream.  Since  the  pressure  is  uniform  across  the  wake,  the  above  immed¬ 
iately  reduces  to 


S  -  Se  »  C  In  (h 
n 

e 


C  In 
P 


iw2 

2  X 


I  + 

2 


1  W2 

2 


Cp  In  (1+x) 


where 

(W-w^)(W  +  w^)/2 
*  "  T  .  i 

i  +  2  -  2 


The  term  x  can  be  considered  small.  Therefore,  the  expansion 


ln(l+x)  “x-YX^+^x^ 


is  approximately  equal  to  x  and  the  expression  for  entropy  becomes: 


(27) 
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Here,  the  approximation  that 


JsCW-Kj^)  -  W 


has  already  been  made. 


Now  following  the  same  steps  as  were  made  earlier  for  the  stagnation 
enthalpy  increment,  it  follows  that: 


Q 

\f 


Mg  6  Sg  -  At  h^  /  (S  -  Sg)  p  w^  dy 


At  h  p 

to  ®  n 


o  O 

W**  /*  /  w  .  pw  , 

J  \  W^'p^W 


M  6S  -  At  h  p  -t —  0 
g  g  m  e  h^ 


Finally,  the  passage-mass-weighted  entropy  increment  is 


M  6S 


8  8  _  JL  JL 

M_  R  Y-1  XC  h  Zirr 
i  n  6 


where  R  is  the  gas  constant  and  y  is  the  ratio  of  specific  heats. 

Equations  (26)  and  (29)  are  expressions  for  the  energy  and  entropy  carried 
by  the  aggregate  wake  particle.  The  effective  streamwise  velocity  at  which  it 
travels  is  now  to  be  defined. 


The  mass-energy-weighted  average  streamwise-velocity-decrement  of  the 
wake,  ug,  is  given  by 

MgAHgUg  -  At  h^  y(W  -  w^)  (H  -  H^)  pw^  dy 
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which,  using  Equation  (27)  and  the  formula  of  Section  3.1.2,  leads  to 


6 

-‘y 

P 


or 


D  G 


13 


D  G 


12 


(30) 


Behind  a  rotor,  the  fluid  In  the  wake  Is  hotter  than  the  free-stream  fluid 
because  more  work  has  been  done  by  accelerating  it  to  a  higher  swirl  velocity. 
Behind  a  stator,  the  argument  Is  rev^'rsed.  The  wake  is  assumed  to  possess  the 
same  stagnation  enthalpy  as  the  free-stream  fluid  and  the  stagnation  rothalpy 
Increment  is  determined.  Writing  again  the  definition  for  total  rothalpy,  for 
any  particle  within  the  wake: 

I“H-Uc®*H-  Uc^  sin  a 

and  for  the  frees tream 
Ig  «  H  -  UC  sin  a 
The  energy  increment  is  then 
I  ••  Ig  ■  U(C  -  Cx)  sin  a 


(31) 

■  At  h  p  UC^  sin  a  0 
m  e 


and  the  mass  integrated  energy  increment  is 


-5 


e'  ^  ~x 
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■  At  h  p  UC  sin 
m  e 
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Finally,  the  passage  mass-weighted  energy  increment  is: 


AI 


M  6l 
_S_8. 


“t 


Ng  UC^  sin  a  0 
Zirr  A  C 

m 


(32) 


Notice,  that  Equation  (32),  for  the  stator  trailing  edge,  is  very  similar  to 
Equation  (26)  which  was  derived  for  the  rotor. 


3.4.3  Radial  Transport  of  Aggregate  Wake  Fluid 

The  radial  movement  of  the  aggregate  wake  fluid  across  the  mean  stream¬ 
lines  is  simply  the  integral  of: 


Dr 


u 

Dx 

w 

g 


(33) 


where 


Dx  ■  differential  distance  in  streamwise  direction 

Wg  •  velocity  in  the  streamwise  direction  of  the  aggregate  wake  fluid 
(see  Figure  3) 

and  the  incremental  radial  velocity  is  reduced  from  the  peak  value  by  the  same 
ratio  as  the  streamwise  aggregate  to  peak  velocity  ratio  [see  Equation  (30)  of 
Section  3.4.2]. 
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(34) 
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Instead  of  integrating  Equation  (33),  a  transformation  is  made  to  the  stream 
function  coordinates,  i)»: 


W  .  Dtp  «  2TTr  X  p  C  cos  4*  Dr 
tt  g  m 


(35) 


where  tg  is  a  dimensionless  stream  function  and  W^-f  Is  the  flow  rate  In  the 
annulus.  Substitution  of  Equation  (34)  into  Equation  (35)  leads  to: 


gte 


Kt  f  ‘ 

TE 


2TTr  X  p  Cjn  u^-g  cos  (p  Dt 


(36) 


with 

Dt  “  Dx/w 

g 


The  values  of  tpg  at  the  trailing  edge  are  taken  as  the  calculation  grid  values 
at  the  trailing  edge  of  the  shedding  blade  row.  At  this  station  tpg  varies  be¬ 
tween  zero  and  unity. 

The  difficulty  that  now  arises  is  that  the  integrated  values  of  i|»g  quickly 
go  out  of  the  zero-to-unity  bounds.  The  formulas  for  the  radial  velocity  G 
take  no  account  of  the  presence  of  the  end  walls,  because  of  the  omission  of 
constraints  imposed  by  a  continuity  equation  and  the  local  change  in  radial 
pressure  gradient  which  must  exist  to  retard  the  radial  wake  movement. 

Instead,  this  problem  is  handled  by  (1)  simply  letting  the  ipg  values  go 
out  of  bounds  and  (2)  defining  a  buffer  zone  adjacent  to  both  end  walls  into 
which  the  wake  fluid  arriving  outside  the  end  walls  is  relocated.  This  pro¬ 
cess  is  illustrated  in  Figure  4.  The  adjustment  is  taken  to  be  quadratic  in 
form,  allowing  a  smooth  transition  of  particle  displacement  at  the  buffer 
boundary.  The  adjusted  values  of  tg  are  denoted  as 


3.4.4  Spanwise  Redistribution 

To  illustrate  the  process  by  which  energy  and  enthalpy  are  transported 
radially,  the  passage-mass-weighted  property  increment,  4P,  is  plotted  versus 
in  Figure  5  at  the  shedding  blade  row  trailing  edge  and  at  some  downstream 
station.  For  simplicity,  aP  will  be  taken  as  any  of  tne  following: 
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Predicted  Movement  of  Wake  Fluid  Aggregates. 


Figures  8  and  9  compare  measurements  of  streamwise  and  radial  wake  veloci¬ 
ties  as  functions  of  downstream  distance,  with  the  values  predicted  by  the  for¬ 
mula  of  Section  3.0.  One  requirement  needed  before  prediction  can  be  made  is 
the  location  of  the  virtual  origin  of  the  decay.  (The  virtual  origin  enters 
the  expressions  as  a  function  of  the  wake  width  and  thus  should  be  the  same  for 
both  streamwise  and  radial  components.)  The  value  selected,  Xg  *  0.12  x 
(chord),  was  chosen  to  be  the  best  compromise  for  this  machine  at  a  C^/Ugi  value 
of  0.85.  At  the  time,  it  was  felt  that  this  value  of  Xg  would  be  changed  as 
experience  was  gained  with  other  test  vehicles,  but  in  practice  this  has  proved 
not  to  be  the  case  and  the  value  of  0.12x  (chord)  is  built  into  the  coding. 

Agreement  between  measurement  and  theory  is  less  favorable  for  the  decay 
of  the  streamwise  deficit  at  a  of  0.75.  The  measured  decay  of  the  deficit 

is  considerably  more  rapid  than  the  prediction  suggesting  that,  perhaps,  for 
this  more  highly  loaded  case  the  assumptions  made  in  the  development  of  the 
theory  are  not  sufficiently  stringent  or  that  the  blade  boundary  layer  was 
separated.  Agreement  of  the  radial  outward  decay  is  good  in  all  cases. 


4.1.2  The  ONERA  Rotor  (Larguier,  References  16  and  17) 

This  rotor  was  used  as  a  calibration  case  by  Adkins  and  Smith  (Refer¬ 
ence  1)  where  the  main  item  of  interest  was  the  peak  radial  velocity.  In  the 
current  investigation,  attention  has  been  drawn  to  the  decay  of  the  streamwise 
velocity  defect  at  three  values  of  radius  ratio. 

These  results  are  shown  in  Figure  12  of  Reference  17.  A  comparison 
between  theory  and  experiment  for  each  radius  is  shown  here  in  Figure  10. 

(The  experimental  values  of  the  defect  were  obtained  as  described  by  Dring  in 
Reference  15.)  Again,  the  value  of  Xg  was  that  used  originally  in  the  anal¬ 
ysis  of  the  Dring  rotor;  it  can  be  seen  that  agreement  between  theory  and 
experiment  is  good  at  each  radius.  One  impressive  feature  of  this  rotor  is 
the  agreement  obtained  between  measurements  made  with  different  forms  of 
instrumentation  (hot  wire,  pressure  probe,  and  laser  velocimeter ) .  For  this 
reason,  it  is  felt  that  matching  this  data  is  a  significant  step  in  the  vali¬ 
dation  of  the  wake  model. 

At  first  sight,  the  peak  radial  velocity  calculated  by  the  new  technique 
at  approximately  midspan  (Figure  11)  appears  to  be  in  close  agreement  with 
that  calculated  by  Adkins  and  Smith.  It  should  be  borne  in  mind,  however, 
that  the  Adkins  and  Smith  peak  radial  velocity  is  calculated  a  short  distance 
(typically  20Z  axial  chord  projection)  downstream  of  the  trailing  edge,  whereas 
the  new  method  calculates  a  trailing  edge  value.  At  the  measurement  station 
(not  shown),  it  is  found  that  the  Adkins  and  Smith  value  is  closest  to  that 
measured  by  a  hot-wire  probe,  whereas  the  new  technique  is  in  best  agreement 
with  rotating  pressure  probe  data. 
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Immersion 


Figure  6.  Drlng  Rotor  -  Comparison  Between  Calculated  and  Measured 
Trailing  Edge  Momentum  Thicknesses. 


4.0  RESULTS 


The  results  obtained  during  the  course  of  this  program  fall  into  two  cat¬ 
egories:  those  concerned  with  the  development  of  the  blade  wake  model  as 
such,  and  those  in  which  the  inclusion  of  the  wake  properties  into  the  mixing 
process  are  compared  with  the  results  obtained  using  the  original  method  of 
Adkins  and  Smith  (Reference  1). 


4.1  WAKE  DEVELOPMENT 


We  are  concerned  with  (1)  the  successful  prediction  of  wake  properties 
(thickness,  streamwise  and  radial  velocity  defects)  at  the  trailing  edge  of 
the  shedding  blade  row,  and  (2)  the  successful  prediction  of  the  variation  of 
these  properties  with  downstream  distance. 

The  test  cases  that  have  been  run  to  verify- the  wake  model  are  discussed 
in  References  15  through  19.  In  general,  each  consists  of  a  low-speed  rotor 
whose  do%mstream  flow  characteristics  have  been  measured  with  stationary  and/or 
rotating  hot  wire  anemometry  and,  in  some  instances  (notably  References  16  and 
17),  with  stationary  and  rotating  pressure  probes  as  well.  Data  are  presented 
at  various  locations  downstream,  thereby  providing  many  checks  on  the  model. 
Questions  have  been  raised  concerning  the  origin  of  the  distance  measurements 
for  References  18  and  19. 


4.1.1  The  Dring  Rotor 

This  rotor  Is  described  in  Reference  15.  It  consists  of  an  isolated 
rotor  in  a  cylindrical  duct.  Predicted  results  have  been  compared  with 
midspan  data  at  four  downstream  locations  under  three  different  flow  condi¬ 
tions  defined  by  C^^/Uq,  the  ratio  of  the  midspan  axial  velocity  to  the  midspan 
wheel  speed  of  the  rotor. 

Figure  6  has  been  included  for  interest.  It  shows  a  comparison  between 
measured  and  calculated  trailing  edge  momentum  thickness  versus  diffusion 
parameters  at  values  of  0.75,  0.85,  and  0.95.  It  can  be  seen  that 

at  Cj^/Uq  ■  0.95  the  calculated  diffusion  parameter  is  somewhau  less  than 
that  measured,  but  at  all  three  conditions  the  momentum  thickness  is  very  close 
to  the  error  band  of  the  measurements. 

Figure  7  compares  the  results  obtained  from  the  different  methods  of  cal¬ 
culating  the  peak  radial  velocity  in  the  wake  -  those  used  by  Adkins  ''-'th 

and  the  new  technique  described  in  Section  3.0.  Again,  these  were  for  _ _ 

Dring  rotor  at  a  Cx/Ugj  of  0.85.  It  can  be  seen  that  the  Adkins  and  Smith 
oiodel  gives  higher  values  than  the  other,  especially  at  the  tip.  At  the  hub, 
the  effect  of  blade  loading  becomes  apparent  with  the  new  technique  showing  an 
increasing  trend,  while  the  Adkins  and  Smith  results  decrease  with  decreasing 
radius.  It  is  felt  that  the  new  method  should  model  the  real  flow  more  closely 
than  the  original  method. 
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which  is  in  agreement  with  Equation  (39.1)  of  the  previous  section. 

Conversely,  if  the  radial  migation  term  Ah',  is  zero  -  either  because  (1) 
the  wake  radial  velocity  is  zero  or  (2)  both  the  energy  increments,  AH,  and 
the  radial  velocities  are  radially  uniform  -  then  the  right-hand  side  is  zero 
and  the  mixing  equation  reverts  to  the  Adkins  and  Smith  form. 

The  formulation  of  the  numerical  solution  to  Equations  (40),  (41),  and 
(42)  is  quite  complex,  partly  because  the  calculation  stations  in  the  through- 
flow  analysis  program  may  not  be  orthogonal.  The  details  of  the  exact  algor¬ 
ithm  are  presented  in  Sections  7.0  and  8.0. 


(40) 
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The_computational  procedure  is  as  follows:  values  of  e,  AH  ,  AX  ,  AS  , 

<0,  and  u  are  computed  by  the  formulas  of  Sections  3.1  to  3.4  of  this  report 
and  Sections  3  and  4  of  Reference  1.  These  calculations  are  uncoupled  from  the 
through-flow  program. 

The  values  of  these  six  parameters  (or  their  equivalents)  at  each  stream¬ 
line-  calculation  station  node  point  are  then  input  to  the  through-flow  program. 
This  program  combines  the  numerical  solution  to  Equations  (40),  (41),  and  (42) 
with  the  continuity  and  radial  equilibrium  equations  to  obtain  the  circumfer¬ 
ential  average  flow  field.  Current  practice  is  to  specify  the  blade  work  dis¬ 
tributions  through  the  specification  of  (the  best  estimate  of  primary  flow 
blade-to-blade  spouting  angles  which  are  in  turn,  modified  by  the  00  secondary 
flow  under/over  turning  angles  mentioned  above.  During  the  convergence  of  the 
through-flow  solution  these  six  parameters  are  held  constant.  After  converg¬ 
ence,  the  six  parameters  are  recomputed  and  the  process  repeated.  The  total 
solution  is  obtained  by  an  iteration  between  the  through-flow  analysis  module 
and  the  secondary  flow/wake  centrifugation/loss  prediction  module. 

It  should  be  pointed  out  that  in  the  through-flow  module.  Equation  (41) 
is  ignored  within  rotors  (that  is,  the  first  station  aft  of  the  leading  edge 
through  the  trailing  edge)  and  Equation  (42)  is  ignored  within  stators.  At 
these  stations  the  swirl  is  obtained  from  the  prescribed  swirl  angle.  Then 
either  H  or  I  is  found  by  the  relation 

0  rCu  -  H  -  X  (43) 

Notice  that  in  the  event  that  e  ■  0,  the  above  equations  reduce  to,  for 
example. 


mf  ..V  ■.*  'w"  •-*  ■  ■ 


where,  as  discussed  in  Section  3.4.1,  the  two  values  of  4p^  and  4Pg  are  equal 
and  opposite  and  therefore  do  not  contribute  to  the  equation.  At  some  down¬ 
stream  station,  however,  the  wake  and  freestream  increments  of  P  have  redis¬ 
tributed  themselves  and,  for  this  same  streamline,  are  denoted  by  4P^  and 
APfs-  So 

P  -  Pte  ♦  AP^  +  APf3 
The  streamwise  change  in  P,  then,  is: 

-5  D  AP  D  AP: 

DP  ^  _ w  fs 

Dm  Dm  Dm 


But  again,  as  indicated  in  Section  3.4.1,  the  change  in  the  freestream  AP  is 
presumed  to  be  much  smaller  than  the  wake  AP  and  is  neglected.  Additionally, 
subscript  w  is  dropped.  The  result: 

DP  -  D  AP' 

Dm  Dm  (39 . 1) 


3.5  NEW  MIXING  SCHEME 

3.5.1  Inclusion  of  Wake  Transport  Terms  in  the  Adkins  and  Smith  Mixing 
Calculation 

Mixing  models  have  now  been  developed  for  complementary  cases: 

1.  Radial  diffusion  of  property  P(r)  with  no  allowances  made  for  cir¬ 
cumferential  variations  of  P  (that  is,  AP  ■  0) 

2.  Radial  migration  of  local  "hot"  wake  fluid  elements  characterized  by 
increment  AP'(r)  referenced  to  a  radially  uniform  free-stream  value, 

Pe  =  F. 

The  general  case  of  radially  varying  P(r)  and  AP  (r)  is  handled  by  "add¬ 
ing"  the  two  procedures  together.  There  is  assumed  to  be  no  significant  mutual 
interaction.  On  this  basis,  the  following  three  equations  are  written: 


Then  by  substitution  of  Equation  (38)  into  Equation  (39)  and  by  using 

d\t  “  2-nr  X  p  C  dn 
m 

re  have,  for  either  side  of  Equation  (39): 
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Hence,  the  mass-weighted  value  of  (P-Pg)  is  held  invariant.  This  is  the 
desired  result. 

The  numerical  method  for  forcing  these  integrals  to  be  invariant  is: 

1.  Fit  a  staggered  parabolic  spline  to  the  curve  of  dPj  versus  4igj  at 
the  shedding  blade  row  trailing  edge,  then  integrate  the  spline  to 
determine  each  individual  area,  SDPj. 

2.  At  downstream  stations,  find  the  staggered  spline  coefficients  which 
yield  the  prescribed  values  of  SDPj  versus  • 

3.  From  the  coefficients  determined  by  (2),  interpolate  for  the  APj 
values  at  the  through-flow  analysis  stream  functions,  j . 

These  AP'  values  represent  the  redistributed  passage-mass-weighted  wake  prop¬ 
erty  increments  for  each  through-flow  streamline  and  station  downstream  of  the 
shedding  row. 

The  change  in  circumferential  average  values  can  be  computed  as  follows: 
before  the  mixing  process  starts,  at  the  blade  t'  >  ing  edge,  the  average 
value  of  P  on  some  given  streamline  is: 

P  -  P_  +  AP  .  +  AP,,. 


or  in  general : 
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where  P  is  stagnation  enthalpy,  stagnation  rothalpy,  entropy,  or  QrCu. 

The  point  of  Figure  5  is  that,  as  the  fluid  moves  downstream  and  the  wake 
moves  radially,  the  individual  areas  under  the  curves,  marked  1,  2,  3,  and  4 
are  held  invariant.  By  assumption,  the  dotted  line  boundaries, 

’^'-•■1/2*  taken  as  halfway  between  the  node  points,  4>v,  ^ 


The  justification  for  this  is  as  follows:  if  we  require  that 
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4.1.3  Penn  State  Compressor  (Reference  18) 

The  Penn  State  University  (PSU)  compressor  consists  of  an  inlet  guide 
vane  (IGV)  row  followed  by  a  rotor  with  a  stator  row  a  significant  distance 
downstream.  The  stage  was  originally  designed  (with  more  traditional  spacing) 
by  Smith  in  the  early  1950* s.  Since  its  migration  to  Penn  State,  a  great  deal 
of  data  have  been  taken  behind  the  rotor  using  both  stationary  and  rotating 
instrumentation.  For  the  purposes  of  calibrating  the  constants  in  the  wake 
model,  it  is  felt  that  data  from  stationary  probes  could  be  misleading  because 
of  the  presence  of  IGV  wakes  which  have  been  convected  through  the  rotor. 
Consequently,  as  far  as  possible,  only  data  from  rotating  instrumentation  have 
been  considered. 

The  results  obtained  are  shown  in  Figures  12  through  15.  Figures  12  and 
13  respectively  show  the  radial  variations  in  peak  spanwise  velocity  and  maxi* 
mum  streamwise  velocity  defect.  Also  shown  are  experimental  data  from  Refer¬ 
ence  18  taken  at  12Z  true  chord  axially  downstream.  Figure  12  shows  that,  in 
this  case,  both  available  methods  for  calculating  peak  radial  velocity  show 
the  same  trend  with  radius  (in  contrast  with,  for  example.  Figure  11). 

Figure  14  shows  the  variation  in  streamwise  velocity  defect  with  down¬ 
stream  distance  at  three  values  of  radius  ratio.  The  hub/tip  ratio  of  the 
machine  is  0.5,  so  these  represent  immersions  of  40.5%,  54.1%,  and  68.4%  of 
annulus  height,  respectively.  It  appears  that  in  all  cases  the  trailing  edge 
defect  is  underestimated.  At  the  outer  two  radii,  the  calculated  and  measured 
defects  are  in  good  agreement  at  streamwise  distances  greater  than  approxi¬ 
mately  20%  chord.  At  the  innermost  radius  (R  =  0.6581),  agreement  is  not  so 
good. 


Figure  15  shows  the  streamwise  decay  of  wake  radial  velocity  at  the  same 
three  radii.  The  calculated  trailing  edge  values  have  been  obtained  using  the 
new  technique,  idiile  the  data  are  taken  from  Reference  18.  It  can  be  seen 
that  the  calculation  seriously  underestimates  the  reported  results.  It  was 
for  this  particular  test  case  that  Adkins  and  Smith  (Reference  1)  reported  the 
greatest  discrepancy  between  measurement  and  calculation.  One  possible  expla¬ 
nation  may  lie  in  the  unconventional  shape  of  the  blades,  that  is,  the  "hooked" 
trailing  edge  of  the  NACA  AlO  meanline  tending  to  encourage  centrifugation. 

One  other  point  raised  by  Adkins  and  Smith  concerns  the  efficiency  at  which 
the  rotor  was  running.  This  was  originally  reported  as  around  82%  and  later 
revised  to  about  86%,  both  of  which  are  low  for  the  condition  of  "good  effi¬ 
ciency"  applied  to  the  calculation  procedure,  and  much  lower  than  measured  by 
Smith  in  the  1950' s. 


4.1.4  Penn  State  Fan  (Reference  19) 


Another  test  vehicle  from  the  PSU  stable,  this  fan  consists  of  12  blades 
of  zero  camber  running  some  distance  downstream  of  a  row  of  support  fins  in  a 


TE  Streamwise  Distance/Chord 


Figure  15.  PSU  Rotor  -  Decay  of  the  Peak 
Radial  Velocity  in  the  Wake. 


cylindrical  duct.  Data  have  been  taken  with  rotating  instrumentation  at  vari¬ 
ous  blade  incidences;  Figures  16  and  17  compare  measurements  and  calculation 
at  a  midspan  incidence  angle  of  10*  and  a  radius  ratio  R  *  0.721  (which  is  the 
midspan  position).  Once  again,  agreement  between  calculation  and  measured 
data  appears  good  for  both  peak  radial  velocity  and  streamwise  velocity 
defect.  Because  of  the  normalization  used  in  the  measured  data,  the  radial 
velocities  shown  in  Figure  17  have  been  plotted  against  normalized  axial 
distance.  Also,  it  is  worth  noting  that  the  two  furthest  downstream  data 
points  in  this  figure  were  taken  with  a  stationary  probe. 

4.1.5  General  Comments  on  Wake  Prediction  and  Decay 

The  value  of  0.12  x  (chord)  upstream  of  the  trailing  edge  for  the  virtual 
origin  of  the  wake  decay  was  chosen  originally  (Reference  11)  as  a  compromise 
between  those  values  that  fitted  the  wake  defect  decay  data  and  peak  radial 
velocity  decay  data  of  the  Dring  rotor  when  it  was  run  at  a  (axial 

velocity /midspan  rotational  speed)  of  0.85.  It  was  anticipated  that  this 
value  would  be  changed  when  other  data  were  available  for  comparison  with  the 
calculated  results,  but  this  appears  to  be  a  reasonable  compromise  for  both 
radial  and  streamwise  decay  of  the  wakes  in  the  low-speed  machines  reported  on 
here.  It  would  be  interesting  to  obtain  data  from  more  realistic  test  vehicles 
in  the  future  in  order  to  verify  this  choice. 

Areas  that  may  require  calibration  in  the  peak  radial  velocity  analysis 
of  Section  3.0  are:  (1)  the  selection  of  the  immersion  of  the  representative 
point  in  the  blade  boundary  layer  (currently  52  of  the  boundary  layer  thickness 
away  from  the  surface)  and  (2)  the  fashion  in  which  the  boundary  layer  form 
factor  changes  from  a  selected  value  (currently  the  flat  plate  value)  at  mid¬ 
chord  to  the  value  given  by  the  Koch  and  Smith  correlations  (Reference  14)  at 
the  trailing  edge  (currently  parabolic).  Again,  for  the  results  presented 
here,  no  changes  have  been  made  to  the  original  scheme. 

With  the  notable  exception  of  the  radial  velocity  in  the  wake  of  the  Penn 
State  rotor  (Section  4.1.3  above),  it  appears  that  the  experimental  results  and 
the  calculations  are  in  excellent  agreement.  With  the  same  exception,  this 
agreement  seems  least  good  for  the  wake  speed  defect  of  the  Dring  rotor  at 
^x/^m  ”  0.75,  followed  by  the  streamwise  decay  of  the  velocity  defect  of 
the  Penn  State  rotor  at  a  radius  ratio  R  ~  0.6581.  These  cases  are  not  con¬ 
sistent,  however  .  For  each  flow  rate  examined  for  the  Dring  rotor,  calcu¬ 

lated  defect  is  higher  than  that  measured;  the  conve  se  is  true  of  the  Penn 
State  data.  It  would  be  tempting  to  suggest  that  the  calculation  of  the  trail¬ 
ing  edge  value  of  the  velocity  defect  needs  re-thinking  were  it  not  for  the 
excellent  agreement  obtained  with  the  ONERA  rotor  and  Penn  State  fan  data 
(References  17  and  19,  respectively). 

Overall,  it  is  felt  that  the  new  approach  to  the  calculation  of  the  trail¬ 
ing  edge  values  of  the  wake  velocity  defects  in  both  streamwise  and  radial 
directions  and  their  decay  with  distance  downstream  represents  a  significant 
improvement  over  the  Adkins  and  Smith  technique. 


4.2  MIXING  CALCULATION  RESULTS 


In  this  section  results  are  presented  for  the  total  mixing  model  as  it 
is  incorporated  in  the  General  Electric  Circumferential  Average  Flow  analysis 
with  Mixing  (CAFMIX)  program.  Comparisons  are  made  with  the  Adkins  and  Smith 
model  and  with  the  case  of  no  spanwise  mixing.  The  vehicle  for  this  presen¬ 
tation  is  the  Air  Force  1500  ft/s.  Transonic,  High-Through-Flow,  Single-Stage 
Axial  Flow  Compressor,  test  data  for  which  are  reported  in  Reference  20.  Fig¬ 
ure  18  illustrates  the  flowpath  and  denotes  the  test  data  point.  Also  shown 
is  the  numerical  streamline/station  grid  used  in  the  calculation. 

A  data  match  analysis  for  this  machine  was  first  performed  by  Simonson 
(Reference  21)  in  1977.  The  ''.urrent  results  extend  this  original  work  to 
include  the  effects  of  the  spanwise  redistribution  of  stagnation  properties 
and  secondary  flow  turning. 

The  principal  new  finding  is  a  more  reasonable  spanwise  distribution  of 
the  rotor  loss  coefficient.  Near  the  hub,  Simonson  found  that  a  very  small 
loss  coefficient  value  of  0.003  was  necessary  to  match  the  experimental  meas¬ 
urements.  With  spanwise  redistribution  of  entropy  by  the  method  of  Section 
3.0,  this  minimum  loss  coefficient  is  a  more  reasonable  0.04. 

With  the  inclusion  of  spanwise  diffusion,  it  is  impossible  to  exactly 
match  the  experimental  distributions  of  total  temperature  and  total  pressure 
reported  in  Reference  20.  This  is  because  the  diffusion  of  these  quantities 
by  Equations  (40),  (41),  (42)  inherently  smooths  the  rotor  discharge  values 
during  their  transport  from  Station  16  to  Station  25.  The  objective,  then'j 
is  to  "match"  the  experimental  data  in  some  smooth-sense.  Figures  19  and  20 
show  how  closely  the  data-match-predicted  total  temperatures  and  pressures 
agree  with  measurements  at  Station  25,  where  the  experimental  results  are 
believed  to  be  reliable. 

Losses  used  in  the  present  effort  are  based  on  the  Koch  and  Smith  (Refer¬ 
ence  7)  and  Smith  (Reference  22)  models.  The  profile  loss  for  both  the  rotor 
and  stator  are  computed  by  the  formulation  of  Koch  and  Smith.  The  end  wall 
losses  are  computed  by  the  method  described  in  the  Appendix  of  Reference  1, 
except  that  the  end  wall  boundary  layer  thickness  in  Equation  (47)  of  Refer¬ 
ence  1  is  reduced  by  using  the  De  Ruyck  and  Hirsch  formulas  (Equations  50,  51, 
and  52b  of  Reference  1).  This  modification  to  the  Smith  repeating-stage  form¬ 
ula  was  made  because  recent  experience  for  "first  stages"  now  suggests  that 
the  end  wall  losses  are  indeed  less  than  for  the  embedded  repeating  stage  and 
because  the  full  repeating  stage  value  of  boundary  layer  thickness  also 
appeared  to  overpredict  the  losses  measured  in  this  application. 

The  rotor  shock  losses  were  initially  based  on  the  Koch  and  Smith  (Refer¬ 
ence  7)  formulation.  However,  it  was  felt  that  this  model  may  be  underpredic¬ 
ting  the  shock  losses.  Consequently,  the  losses  for  the  outer  portion  of  the 
rotor  have  been  increased,  as  shown  in  Figure  21,  to  yield  a  better  match  with 
the  measurements.  Figures  21  and  22  compare  the  present  loss  coefficient  dis¬ 
tribution  (labeled  by  W.K)  with  those  found  by  Simonson,  which  included  no 
spanwise  mixing. 


Streamlines  and  Calculation  Stations  for  the  Air  Force  High  Through  Flow  Stage 


with  No  Mixing  (Simonson). 


Also  shown  is  the  modified  loss  model  of  Koch  &  Smith  and  Adkins 
(labeled  LOSS.MDL). 


Stator  Logs  Coefficient  Used  in  the  Present  Study  (W.K)  and  in  the 
No-Mixing  Data  Analysis  (Simonson). 


The  stator  losses  used  in  the  present  method  (labeled  W.K  in  Figure  22) 
are  somewhat  less  than  those  ascertained  by  Simonson  from  examination  of  the 
stator  wake  data,  but  the  shape  of  the  curve  is  very  similar. 

The  mixing  phenomenon  upsets  the  normal  method  for  emprically  determining 
losses.  Usually,  the  loss  coefficient  is  determined  from  the  measured  total 
pressure  just  in  front  of  and  behind  the  blade  row  at  immersions  with  the  same 
stream  function.  The  loss  coefficient  computed  by  CAFMIX  in  this  way  leads  to 
the  pseudovalues  shown  in  Figure  23. 

Figure  24  compares  the  data  match  values  of  the  (primary  flow)  rotor  dis¬ 
charge  swirl  angles  with  the  design  angles  reported  in  Reference  20.  Also 
shown  are  two  values  calculated  by  a  coupled-boundary  layer/inviscid  blade- 
to-blade  solution.  These  latter  values  are  in  rather  good  agreement  with  the 
present  analysis.  The  predicted  full  three-dimensional  circumferential  flow 
angles  (primary  and  secondary)  are  illustrated  in  Figure  25.  There  is  little 
secondary  flow  under/overturning  except  at  the  end  walls. 

The  remaining  curves  in  this  section  illustrate  the  intermediate  values 
of  various  parameters  of  the  wake  centrifugation  model  as  applied  to  the  Air 
Force  stage;  they  compare  these  parameters,  where  appropriate,  with  the  Adkins 
and  Smith  model.  Figure  26  shows  the  predicted  rotor  wake  width  as  defined  in 
Figure  27.  Note  that  the  two  models  use  very  different  wake  widths.  Adkins 
and  Smith  base  their  wake  width  on  total  losses  (including  shock  losses  but 
not  end  wall  losses);  they  back-calculate  a  momentum  thickness  and  wake  thick¬ 
ness  as  described  in  Section  4.2  of  Reference  1.  In  the  current  model, 
denoted  as  Whitfield  and  Keith,  the  wake  width  is  consistent  with  the  profile 
boundary-layer-growth  correlations  used  by  Koch  and  Smith  to  evaluate  profile 
loss.  The  boundary  layer  is  then  thickened  to  account  for  one-third  of  the 
profile  shock  loss. 

The  model  for  the  radial  velocity  in  the  wake  has  been  discussed  in  some 
detail  in  Section  3.1.1.  Results  for  several  machines  were  presented  in  Sec¬ 
tion  4.1.  The  improved  model,  as  shown  in  Figure  28,  yields  higher  values  for 
the  rotor  wake  radial  velocity  than  does  the  Adkins  and  Smith  formulation. 

The  radial  velocities  are  small  near  the  hub,  but  become  quite  large  (approx¬ 
imately  70Z  of  the  meridional  through-flow  velocity)  near  the  tip.  Conversely, 
the  stator  wake  radial  velocities  (not  shown)  are  small  across  the  span  (about 
negative  12  of  the  meridional  velocity). 

The  radial  secondary  flow  velocities  behind  the  rotor  are  shown  in  Fig¬ 
ure  29.  They  are  identical  for  the  two  models. 

The  wake  width,  wake  radial  velocity,  and  radial  secondary  velocity  con¬ 
stitute  the  elements  of  the  blade-to-blade  trailing  edge  spanwise  velocity 
field  illustrated  in  Figure  27.  This  curve  is  integrated,  and  the  mixing 
coefficient  is  calculated  by  the  formulation  presented  in  Section  4.0  of 
Reference  1.  The  resulting  mixing  coefficients  are  presented  in  Figures 
30  and  31.  The  Whitfield  and  Keith  coefficients  (Figure  30)  are  initially 
higher  at  the  rotor  trailing  edge)  than  the  Adkins  and  Smith  values, 
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Spanwise  and  Streamwise  Variation  of  Stagnation  Temperature  Increment 
Generated  by  the  Rotor  Wake,  Whitfield  and  Keith. 


Radial  Velocity,  G,  with  Downstream  Distance,  Whitfield  and  Keith 


but  then  fall  off.  This  is  due  to  the  use  of  the  local  decayed  value  of 
radial  wake  velocity,  as  predicted  by  the  formulation  of  Section  3.2, 
rather  than  the  invariant  (trailing  edge)  velocities  used  by  Adkins  and  Smith. 

The  decay  of  the  wake  radial  velocities  for  three  representative  stream¬ 
lines  is  shown  in  Figure  32.  In  Figure  33,  trajectories  of  aggregate  wake 
fluid  particles,  moving  outward  by  velocities  related  to  those  shown  in  Fig¬ 
ure  32,  are  plotted  against  a  background  of  mean  flow  streamlines.  Notice  the 
large  outward  movement  of  the  wake  fluid.  (These  effects  are  computed  by  Equa¬ 
tion  (20)  of  Section  3.2.2  and  Equation  (36)  of  Section  3.4.3,  respectively.) 

In  the  midspan  region  of  the  annulus,  the  radial  distance  between  the 
wake  particles  is  increased  as  the  wake  moves  downstream.  This  thins  the 
wake  and  causes  the  mass-averaged  AT'f  (Figure  34)  and  AS'  (Figure  35) 
to  decrease.  Conversely,  near  the  casing  the  wake  fluid  tends  to  pile  up 
and  AT.j^  and  AS'  increase.  However,  the  hub-to-tip  integrated  levels  of  both 
ATj  and  AS'  remain  constant.  The  effects  shown  in  Figures  33,  34,  and  35  are 
related  to  the  rotor  wake.  Interaction  of  the  rotor  wake  with  the  stator 
blade  flow  field  is  ignored,  as  is  impingement  of  the  wake  on  the  pressure 
side  of  the  stator  blade. 

The  radial  migration  of  the  stator  wake  is  small.  Consequently,  the 
plot  of  AS'  for  the  sequence  of  stations  behind  the  stator  trailing  edge  all 
fall  one  upon  the  other  (Figure  36).  The  methods  for  predicting  AT^  and  AS' 
are  given  in  Section  3.4.3. 

Finally,  for  reference,  the  source  terms  in  the  spanwise  mixing  equations 
(Equations  40  and  42}  are  plotted  in  Figures  37  and  38.  Again,  it  may  be  seen 
that  these  source  terms  must  integrate  (with  respect  to  mean  flow  stream  func¬ 
tion)  to  zero. 

An  overall  comparison  between  the  present  study  and  Adkins  and  Smith  is 
provided  in  Figures  39  and  40.  The  two  predictions  use  the  same  loss  coeffi¬ 
cients  and  rotor  and  stator  (primary  flow)  exit  angles.  The  differences  shown 
indicate  the  differences  that  result  from  (1)  the  revised  mixing  coefficient 
and  (2)  the  inclusion  of  "hot"  wake  radial  migration  effects.  The  differences 
are  not  large,  but  they  are  significant  for  this  case. 

This  significance  is  further  illustrated  in  Figure  41.  This  is  a  plot  of 
the  entropy  distribution  at  a  number  of  stations  within  the  Air  Force  compres¬ 
sor.  Both  the  jump  in  entropy  across  each  blade  row  and  the  redistribution  of 
entropy  at  subsequent  downstream  stations  can  be  observed.  The  final  distri¬ 
bution  at  Station  25  is  shown  to  be  in  reasonable  agreement  with  the  experi¬ 
mental  data.  It  is  through  the  wake  centrifugation  process  that  the  rotor 
losses  (near  the  80Z  immersion  point)  are  "thinned"  and  the  combined  rotor  and 
stator  losses  then  agree  with  the  Plane  25  measurements. 

The  calculations  reported  in  this  section  have  been  carried  out  with  an 
invariant  radial  equilibrium  concept.  The  radial  variation  of  streamline  curv¬ 
ature  and  flow  angle,  together  with  certain  streamwise  derivatives,  was  deter¬ 
mined  by  Simonson  (Reference  21).  His  calculations  included  six  stations  in 
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Wake  Radial  Velocity/Meridional  Velocity 


Figure  27.  Construction  of  Spanwlse  Velocity  Field  From  Secondary  Flow  and 
Wake  Contributions. 


Width/Staggered  Spacing 


Match  (Primary  Flow)  Rotor  Spouting  Angles  Compared  to  the  Design  Flow 


Stator  Pseudo  Loss  Coefficient 


Whitfield  and  Keith. 


Whitfield  and  Keith. 


Figure  39.  Comparison  Between  Adkins  and  Smith/  Whitfield  and  Keith  Predicted 
Total  Temperatures  at  Station  25  for  the  Same  Rotor/Stator  Blade 
Exit  Primary  Flow  Angles  and  Loss  Coefficients. 


Figure  41.  Variation  of  Entropy  from  Rotor  Trailing  Edge  to  Measurements  Station  25 
as  Determined  by  the  Solution  of  Equation  42  of  Section  3.5.1. 


the  rotor  and  six  stations  in  the  stator.  Subsequent  calculations  with  the 
Adkins  and  Smith  and  Whitfield  and  Keith  models  have  been  performed  with  these 
through-the-blade-row  stations  removed.  However,  for  consistency  the  latter 
calculations  have  both  used  the  (presumed  more  accurate)  Simonson  streamline 
curvatures. 


5.0  CONCLUSIONS 


The  peak  wake  radial  velocity  calculation  of  Adkins  and  Smith  has  been 
revised  to  include  the  effects  of  blade  loading.  Universal  velocity  pro¬ 
file  shapes  have  been  assumed  for  both  the  radial  and  streamwise  wake  pro¬ 
files  and  these,  together  with  simplified  momentum  equations  in  the 
streamwise  and  radial  directions  and  the  knowledge  of  the  trailing  edge 
form  factor  have  been  used  to  obtain  the  following: 

•  Maximum  streamwise  deficit  at  the  trailing  edge  and  its  decay 
with  distance  downstream 

•  Decay  of  the  peak  radial  velocity  increment/decrement  with 
distance  downstream 

•  Trailing  edge  wake  thickness  and  its  increase  with  distance 
downstream 

•  Trailing  edge  momentum  and  displacement  thicknesses. 

Wakes  have  been  considered  as  streams  of  energy-carrying  particles  which 
travel  at  mass-energy-averaged  velocities  in  the  radial  and  streamwise  . 
directions.  At  calculation  stations  downstream  of  the  shedding  trailing 
edge,  the  properties  carried  by  the  particles  (total  enthalpy  behind  a 
rotor,  total  rothalpy  behind  a  stator,  and  entropy  everywhere)  are  redis¬ 
tributed  on  the  circumferential  average  streamlines  for  use  in  a  revised 
formulation  of  the  mixing  equation  as  a  part  of  the  circumferential  aver¬ 
age  flow  determination  solution  iterative  scheme. 

Comparisons  between  the  theory  and  available  experimental  data  show  good 
agreement  between  predicted  values  of  both  streamwise  and  radial  peak 
velocity  deficits  in  most  cases.  If  the  shedding  blade  boundary  layer  is 
separated,  the  wake  appears  to  decay  significantly  more  rapidly  than  pre¬ 
dicted. 

The  new  mixing  model  has  been  applied  to  the  Air  Force  High-Through-Flow 
Axial  Flow  stage  where  it  is  found  that  the  apparent  super-high  efficiency 
near  the  hub  of  this  machine  is  due  in  part  to  the  centrifugation  of  wake 
fluid  to  the  outer  region  of  the  annulus. 


APPENDIX  A:  STAGGERED  PARABOLIC  SPLINE  FORMULATION 


In  this  appendix  the  equations  for  fitting  a  staggered  parabolic  spline 
curve  through  a  given  set  of  points  are  fomulated.  This  particular  curve  fit 
is  convenient  because  (1)  either  the  points  through  tdiich  the  curve  must  pass 
can  be  prescribed  or  (2)  the  local  integral  of  the  curve  (between  knots)  can 
be  prescribed.  In  Section  A.l  the  recursion  relation  for  a  spline  passing 
through  a  set  of  points,  u^,  is  developed.  In  Section  A. 2  a  formula  for  the 
integral  across  a  spline  segment  is  presented.  Finally,  in  Section  A. 3  the 
recursion  formula  in  terms  of  known  integral  values,  rather  the  u  values,  is 
found. 


A.l  RECURSION  RELATION  FOR  PRESCRIBED  POINT  VALUES 

Using  the  nomenclature  illustrated  in  Figure  42,  the  parabolic  (spline) 
relations  are  written  for  the  i^^  interval  defined  by  h^. 

^i  "  “i  ""i  *  *  T  **1  (Al) 


S  .  *  c(^  +  di  X 


mi  «  ci  -  di  fi 


“i+l  *  ci  +  di  gi+i 

^ere  ci  and  are  unknown  coefficients,  S  is  the  ordinate  of  the 
spline  and  mi  is  the  slope  at  the  knot. 

Subtract  (3)  from  (4) 

mi+i  -  mi  =  di  (fi  +  gi^j) 

Then 


°'i+l  -  ”>i 
^i  ®i+l 


'i-H  "*i 


^  fSi+l  “i  +  “i+iJ 
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8i 


ETlX 


A.  2  REPLACEMENT  OF  POINT  VALUE  BY  INTEGRAL  VALUE 
Over  the  segment  hi, 


*1+1  ii+i/ 

j  ^1“’*  ■  /  <“1 

-^1  -h  * 


‘*i+i"i  +  ‘i^i+P  h 


"i  <*1+1  *  'i>  i  <*i+ri  *i"i+i)  <»  1+1 

2hi 


<"i+i  -  "i>  <*'i+i  - 

6hi 


+  1  (m^^i  -  mj 


^1  ■  “l"!  *  6^  *'l+l  -  2 


i+1  2  2 


2l  ■  '■l  '“l  +  ?  <"l+l  -  ”l>  '■ll  I  *\+l  "l  -  i  f2^"l+l 


A. 3  MATRIX  EQUATION  IN  TERMS  OF  I  RATHER  THAN  u 
Rewrite  (A12)  in  form  to  be  added  to  (AlO). 


®j  o  o  ^<+1  2  2 

+3ir  t^S  i+1  -  -  f  +  3h~  i  ■*■  ®  i+1^ 

1  J-  (A13; 


21, 


m 


-2u..  ,  + 


1-1 


i-1  3h 


i-1 


cVi  +  2fj.i*l  +  £\.,l 


+  fhTTj  1-1  ■  ^^1-1  *1  ■  ®  i' 


21 


i-1 


1-1 


Add  (AlO),  (A13),  and  (A14) 


2  2  ™i-l 

fSl  +  +  f 


i-1 


■*■  1-1  ■  ^^1-1*1  ■  *1  ■*'  ^®l'’l-l  ■'■  ^^l-l®l'  3hj^_j 


■£(2g\„  -  2£^g,^j  -  f2^  *  3f,2  1  1  j 

3hi)  1 


+  [3f^  -  2£\  3-  3  ^  .  2  ^  I 

3n^  1  i-i 


Simplifying  - 

1  ^i-iVi  I  ^^1-1  +  ^1^ 


^1  ^i-i 

,  «  2  f— i _ i-i-1 

3.1+1  ^ 


m. 


or 


hi-i  ni^.i  +  2  (h^_^  +  h^)  ^ 
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(A14) 


APPENDIX  B:  THE  DIFFUSION  EQUATION  FOR  CIRCUMFERENTIAL  AVERAGE  MIXING 


In  this  section  the  diffusion  equation  for  the  mixing  of  stagnation 
enthalpy,  rothalpy,  and  entropy  is  formulated.  Unfortunately,  this  equation 
is  complicated  by  the  requirement  that  the  computational  grid  used  for  through- 
flow  analysis  is  generally  a  nonorthogonal  one. 


B.l  SECONDARY  FLOW  TRANSPORT 

The  secondary  flow  transport  of  any  circumferentially  averaged  property, 
P,  is  assumed  to  be  perpendicular  to  the  streamlines  in  the  meridional  plane 
and  is  given  by 


(Bl) 


where 


n  •  distance  normal  to  streamlines 
p  ■  density 

m  ■  distance  in  the  meridional  direction 
t  ■  coordinate  in  the  blade-to-blade  direction 
a  *  blade  spacing 

Vjj  »  meridional  component  of  circumferential  average  velocity 

Vg  ■  spanwise  component  of  secondary  and  wake  flow  (in  the  n-direction) 

P  >  circumferential  mass-weighted  average  property.  (In  this  section  the 
overbars  on  P  are  omitted.) 

In  terms  of  the  nonorthogonal  coordinates,  q  and  m,  shown  in  Figure  43, 
the  normal  derivative  is  given  by: 


3n 


1 

cos  ((J>  -  o) 


dq 


-  tan  (4)  -  o) 


3m 


(B3) 
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The  second  term  on  the  right  of  (B3)  is  generally  small  since  (^-o)  is  small, 
and,  if  included,  it  would  make  the  parabolic  marching  solution  procedure 
impossible.  Therefore,  it  is  dropped.  So: 


1  3P 

3n  cos  (<|)  -  a)  3q  (B4) 


B.2  DIFFERENTIAL  FORM  OF  THE  DIFFUSION  EQUATION 

Reference  is  now  made  to  Figure  44.  Considering  the  rectangular  torus 
bounded  by  the  faces  labeled  as  i,  j,  i-1,  j-1,  it  is  required  that  the  net 
flux  of  P  through  all  of  these  faces  be  equal  to  the  source  of  P  integrated 
over  the  volume.  As  noted  in  the  previous  section,  the  flux  normal  to  the 
streamlines  is 


3P  £  3P 

^  3n  cos  (4)  -  a)  3q 


(B5) 


The  transport  of  P  parallel  to  the  streamlines,  carried  by  the  fluid,  is 


pV  P 
m 


(B6) 


times  the  area  normal  to  the  streamlines. 
The  flux  balance  then  is: 


[2iTrXh2AC2  ^os  ((j)  -  o)pV^  P]^  -  [2TTrXh2A52  ^^i-1 


+  [2iTrXh2A52  (4)  -  o)  “  [2iTrX  h2A^2  (4>  -  a)  e  PnU-l 


[2iTrXhj^AC^  e  P^]j  +  [2irrAhj^ACj^  e  P^I 


(B7) 


I  I  2TrrXhj^d5j^  cos  (4>  -  o)  h2  d?2 

M-l  A-1 


where  Sc  is  a  source  term  per  unit  volume. 
To  simplify  this  somewhat,  substitute 

f-  *  2irrXho  cos  (d)  -  o)pV 
z  z  m 


(B8a) 


2TrrXhj^  e 
~  cos  (4)  -  a) 


=  27rrXh2  tan  (4)  -  a)  e 


into  (B7).  Using  (BS),  the  result  is: 


(B8b) 


(B8c) 


Plj  -  [fjiCj  Pl^.^  -f  [gjijj  -  [gjSCj 

-  [8iA5i  P,1j  +  lgi45i  P,lj.i  -  RHS 


where  RHS  signifies  the  right-hand'-side  term  of  Equation  (B7).  Since  ^^2  ia 
the  same  at  faces  i  and  i-1  and  ACi  is  the  same  at  faces  j  and  j**l,  these 
values  can  be  divided  out.  The  first  line  of  (B9)  will  become 


[f2P  +  g2Pq]i  -  [f2P  +  82*’qU-l 


which  in  differential  form  will  be 


ftjp  +  gjp  1 


By  similar  arguments  for  the  2nd  line  of  (B9) ,  the  following  differential 
equation  is  written: 


[f2^  +  82Pq]  -  =  2-nrX  cos  (4)  -  0)hj^h2S^ 


(BIO) 


B.3  TRANSFORMATION  TO  STREAM  FUNCTION  COORDINATE 


A  somewhat  more  convenient  form  of  the  differential  equation  is  obtained 
by  replacing  the  ^2  coordinate  with  the  stream  function  coordinate,  <|).  The 
variable  f2,  defined  by  Equation  (B8a)  relates  d(2  with  dt(). 


dij;  =  2TrrXh2  cos  (<{>  -  a)  PV^  dC2 


^2  ^^2 


Also  notice  that  if  the  continuity  equation  is  satisfied,  which  it  is,  then  f2 
is  constant  in  the  ^-direction.  With  this  in  mind.  Equation  (BIO)  can  be 
divided  by  f2  with  the  following  result: 


h,S 

1  c 


(BID 


where 


27rrA  e 
cos  ((J)  -  a) 


(B12a) 


8s 


tan  (4)  -  o)  e  M  ^2 

cos  (4>  -  a)  PV  f^ 
m  z 


(B12b) 


B.4  DIFFUSION  EQUATIOM  SOURCE  TERM 


The  source  term  in  Equation  (Bll)  is  determined  indirectly.  As  discussed 
in  Section  3.5.1,  wake  centrifugation  effects,  in  the  absence  of  a  spanwise 
variation  of  P  (that  is,  ==  Pq  =  0),  led  to  a  circumferential  average  value 
of  P  equal  to  AP'.  Therefore,  the  differential  equation  must  have  the  form: 


[P 


Vql 


SAP* 

3^1 


(B13) 


so  that  when  Pq  ■  0,  it  follows  that 

3P  _  3AP* 

3^1  ■  3Ci 


Therefore,  the  source  term  by  deducton  is 


S 

c 


p\  3AP' 
hi  3Ci 


Equation  (B13), 


together  with  the  definition 


P 

q 


(B14) 


is  the  form  of  the  diffusion  equation  to  be  solved  in  the  through-flow  analysis 
program.  This  equation  is  linear  for  prescribed  values  of  g  and  g^,  both 
defined  by  Equation  (B12). 
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APPENDIX  C:  NUMERICAL  SOLUTION  TO  THE  DIFFUSION  EQUATION  IN 
NONORTHOGONAL  COORDINATES 


No  claim  is  made  that  the  mathematical  model  described  in  this  report  is 
a  highly  accurate  one.  In  fact,  the  model  is  inherently  an  approximate  one. 
Nevertheless,  we  chose  a  good  quality  numerical  scheme  for  integretating  the 
diffusion  equation,  derived  in  the  previous  section,  because  the  hub-to-tip 
streamline  grid  spacings  may  be  very  nonuniformly  distributed.  In  fact,  a 
high  concentration  of  grid  lines  near  the  hub  and  tip  and  a  sparse  midspan 
distribution  is  preferred.  To  prevent  a  solution  that  is  largely  dependent 
on  the  number  of  spanwise  grid  points,  a  second  order  accurate  algorithm  is 
a  requirement.  Another  requirement  is  global  conservation  of  the  property  P. 

To  meet  these  two  requirements,  a  general  spline  scheme  Is  formulated 
wherein  an  Integral  form,  as  well  as  the  differential  form  of  the  conservation 
equation.  Is  employed.  The  basic  differential  equation,  as  developed  In  the 
previous  section.  Is 


tP  +  SsPql 


3(|; 


[ghi  P^] 


3AP* 

'^1 


(Cl) 


where 


ZTTrX  £ 
cos  (4)  -  a) 


(C2a) 


tan  (4)  -  g)  e 

cos  (4>  -  a)  pVn, 


(C2b) 


and  for  later  reference: 


f  *  ZirrX  cos  (4)  -  a)  P'^m 


(C2c) 
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Equation  (Cl)  is  first  discretized  in  the 

[P  +  S^P,!!  -  [P  +  8,P,li.i  -  e[|;  (8h,P,)]j 


-  (8j_j  hjj  Pq^_j)  '  *Pi  -  ^Pf-l 


For  convenience,  we  have  chosen  AC  *  1*  Note  that  the  distance  between  adja¬ 
cent  grid  points  is 

“i  ■  ®i-l  “  ^li  (C4) 


Because  the  integration  is  between  i-1  and  i,  where  (C4)  applies,  it  is  appro¬ 
priate  to  use  h^i  where  h]^  appears.  Hence,  h]^  £-]^  does  not  appear  in  (C3). 
Also,  the  parameter,  0,  is  the  Crank-Nicholson  weighting  parameter.  Ideally, 
e  is  one  half.  However,  a  value  slightly  larger  than  one  half  is  generally 
chosen  for  stability. 

Rearrangement  of  (C3),  to  put  the  "knowns"  on  the  right-hand  side,  gives: 


®  [aip  ^®^1  ^q^]i 


=  8s  ^ii-1  -  ®>  i  pqi_i)  Ap;  - 


This  equation  is  now  integrated  with  respect  to  ^  using  the  trapezoidal  rule: 

j[^  ^8^1  ^q^]j-l  [lip  ^q^]jj  ~  *  ^8^3^  Pq)^  -  (ghj^  ^q^j-1  (C6) 


The  result  is: 


(P  +  s, 

[P  +  gs 

**•1 

1  Vi.j 

-1  j 

= 

[P  +  g^ 

Phi  •  + 
q  i-1.3 

^  ‘'■’’‘-i.J-i  (C„ 

+ 

AP'  .  + 
i,J 

ap!  .  , 
i.J-l 

-  AP'i  .  - 

i-1.3 

■  ^^i-l,j-ll 

+ 

(1  -  6) 

h-,  . 

11.3 

■  ^li,j-l  ^2^q^l-l,j-l 

Equation  (C7)  provides  one  relationship  between  P  and  Pq.  One  additional  rela¬ 
tion  is  needed  to  complete  the  system.  This  relation  is: 


-  '■j-i  -  ^  ^ 


1  = 
-ij  12 


Equation  (C8),  with  1^1,  is  based  on  a  quartic  polynomial  for  P  between  points 
j-1  and  j.  The  quartic  is  reduced  to  a  parabola  by  setting  X“0.  The  second 
derivative,  Pqq,  is  related  to  tetms  in  the  differential  equation  (C5)  in  the 
following  way: 


S*'!  ^  ^  <**'1  ^  ‘S"!' 


qq  gh 


^  [f  i  <*'’1  ^  ^  <8hi)] 


Substitution  of  (C9)  into  (C8)  yields: 


9q 


dip 


(gh. 


^^J-1 


-  P, 


.j-i  i  <*»i>j-i]  I  -  “ 


(CIO) 


Thus  we  have  three  equations,  (C5) ,  (C7),  and  (CIO),  for  the  system 
matrix  equation.  The  three  unknowns  at  each  mode  are  P^,  P,  and 


To  conserve  property  P,  the  boundary  conditions  are: 

Pq  “  0  for  j  *  1  and  j  ■  NJ 

where  j  *  1  and  NJ  corresponds  to  the  casing  and  hub  boundaries. 

The  parameter  X  appears  in  Equation  (C9).  If  X*0,  this  scheme  is  known 
as  the  Keller-Box  method.  Unfortunately,  for  small  values  of  e  (that  is, 
small  secondary  flow)  the  Keller-Box  method  produces  wiggles.  The  wiggles 
disappear  where  A~l,  in  which  case  the  method  is  a  cubic  spline  method.  Inter¬ 
mediate  values  of  X  are  also  permitted. 
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